Limiting temperature from a parton gas with power-law tailed distribution 
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We combine Tsallis distributed massless partons to an effective thermaf prehadron spectrum by 
fofding. A limiting temperature and a mass spectrum combined of three exponentials emerge by 
this procedure. Meson and baryon resonance spectra have different polynomial prefactors. 
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The idea to treat the bulk of mcsonic and bary- 
onic resonances as a statistical system stems from Rolf 
HagedornQ. The mass spectrum is exponential, multi- 
plied with an originally negative power of the mass, which 
fact gives rise to a maximal, limiting temperature, Th, 
when heating a hadron gas. The exponential factor with 
the main part of the hadron energy, i.e. the rest mass, 
is the most essential feature in this assumption, power 
factors and depending on them the quantitative fit to 
the Hagedorn temperature, Th, vary. A recent compi- 
lation of hadron resonances by Broniowski, Florkowski 
and Glozman shows that this idea works well beyond the 
data used for original fits, although baryons and mesons 
seem to follow a slightly different line|2j. Explanations 
for the origin of this exponential mass spectrum date back 
to the MIT bag model, where this behavior was demon- 
strated by KapustaQ. Bugaev and others point out that 
such a system is a perfect thermostat forcing the same 
temperature to any finite system in thermal contact 0. 
The picture of Hagedorn resonances also fit well to lattice 
QCD results on the equation of stated and is under con- 
sideration in recent microscopic models of quark matter 
rehadronization in heavy- ion collisions 0. 

In this paper we demonstrate how another approach 
may support the occurence of an essentially exponential 
mass spectrum of hadrons. We consider massless thermal 
partons, but with a generalized equilibrium distribution 
with power-law tail asymptotics. This distribution will 
then be folded to mcsonic or baryonic one-particle distri- 
butions of energy, or - at zero rapidity - transverse mass 
and transverse momenta, respectively. Such stationary 
distributions are conform with the basic laws of statis- 
tical physics and may be considered as a description of 
the intermediate Opt ~ 1 — 4 GeV) part of the observed 
hadron spectra0,[a]. This approach delivers qualitatively 
interesting results, in particular a characteristic differ- 
ence between the mesonic and baryonic mass spectrum. 

There are numerous occurrences of power-law tailed 
statistical distributions in nature. In particular hadron 
transverse spectra stemming from elementary particle or 
heavy ion collisions can be well fitted at mid rapidity by a 



formula reflecting 



m T -scaling H i, [3 13 13 13 13 : 

p T )~(l + m T /E c )- v . (1) 



Interpreting these spectra in terms of the single particle 
energy, one considers E — tut — \/y\ + for a rel- 
ativistic particle with mass m. This formula describes 
a Tsallis distribution, which was conjectured earlier by 
using axiomatic thermodynamical arguments ^(|. This 
differs from the traditional interpretation of such spec- 
tral tails in particle physics, when these are treated as 
non-equilibrium phenomena. The very high-px part, ex- 
pected to stem from jet fragmentation, may still allow 
for a statistical interpretation below px = 6 — 8 GeV. 

Distributions with a power-law like tail are encoun- 
tered in several different statistical 

models[l[3I3I3 

13 HHH They are investigated as generic distri- 
butions in the non-extensive thermodynamics [3 13113' 
based on a generalization of the Boltzmann-entropy, en- 
countered first in informatics 

problems [3 El- Without 
being able to exclude the non-equilibrium interpretation 
of the power-law tail in particle physics, we explore here 
some consequences of a stationary state with Tsallis dis- 
tributed extreme relativistic particles (massless partons). 

We consider a massless parton gas with binary colli- 
sions obeying a general, non-extensive energy composi- 
tion rule[29j, -E12 = h(E\, £2). Whenever this rule is as- 
sociative, the one-particle energy can be mapped onto an 
additive quasi-energyX (£"12) = X(Ei)+X(E 2 ), with the 
help of a strict monotonic function, obeying A(0) = 0. 
For the Tsallis stationary distribution one uses 



X(E) = - ]n(l + aE), 



(2) 



with an energy scale a = 1/E C related to the micro- 
scopical dynamics. This gives rise to the following non- 
extensive composition rule: 



h(x, y) = x + y + axy. 



(3) 



The stationary one-particle distribution under these con- 
ditions becomes 



f{Ei) = 



1 



ZiW) 



expC-X^/T), 



(4) 
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with a temperature T determined by the conserved to- 
tal (quasi) energy X(Etot) and particle number. The 



2 



onc-particlc partition function here is given by the phase 
space integral Zi((3) — J dTjf(Ej). 

In this paper we are interested in the energy distri- 
bution of large subsystems, i.e. in the microcanonical 
distribution of the energy of N particles each following 
a Tsallis distribution given by eq.Q). The general for- 
mula for the TV-particle energy distribution, assuming a 
factorization of the one-particle distributions is given by 



)• 



r N N 

F N (E) = AEX'(E) / l[dT l S(X(E)-Y,X(E l )) f(E, 

J i=l i=l 

(5) 

Here X'(E) stands for the derivative of the strict mono- 
tonic mapping function X(E), AE is the width of the 
iV-particle energy shell and the dTi integration measures 
refer to the one-particle phase space factors. The factor- 
ization is usually a good approximation for values of N 
being still negligible besides the total number of parti- 
cles. A check of this formula for N = 1 expresses the 
one-particle energy distribution as being proportional to 
the one-particle phase space distribution: 



F^E) = ^E 2 AEf(E). 



(0) 



Since F N (E) contains an X'(E) factor, F N (E)dE = 
gx(X)dX relates it to the distribution of the N— particle 
quasi energy, X = X(E). As it was shown in Ref.p^. 
in certain kinetic models leading to a stationary state of 
the non-extensive thermodynamics, X(E) = ^2nX(Ej) 
is the conserved quantity and the quasi-energy can be 
regarded as the physical energy of composite N— particle 
systems. 

Using the Fourier-expansion of the constraint on the 
sum of the quasi-energies, Xi = X(Ei), the expression 
(JSJ) factorizes: 



g N {X) = AE 



ds 
2^ 



isX 



N 

n 



dT 3 f{E 3 )e 



-isX(Ej) 



(7) 

Utilizing now the equilibrium one-particle quasi-energy 
distribution, Q we obtain 



-isX(Ej) 



Z^P + is) 



(8) 



The TV-particle quasi-energy distribution we arc seeking 
for is then given in a form normalized to one in an energy 
shell of width AE as 



sA 



ZiQg + ta) 

MP) 



N 



(9) 



Such integrals may be evaluated in a saddle point approx- 
imation, which is good for large N as long as no singular- 
ity has been encountered in the expansion of In Z\((3+is). 
The result is a Gaussian distribution in X{E). Irrespec- 
tive to this approximation, as long as the factorization 



assumption is valid, the exact expectation value is given 
by 



and the square width by 



if 



5X(EY = AT—lnZrC/J) 



(10) 



(11) 



In case of the Tsallis distribution using J5J for massless 
particles one obtains 

v ' fe=i 

with spatial volume Vd in d dimensions and with the fac- 
tors Bk = P — ka. The expectation value of the quasi- 
energy per particle becomes 



ei = 



( X(E)} 
N 



(13) 



k=l 



while the unit square width contribution is given by 

(14) 



SX 2 x ^ 

k=l 



All these expressions (I12I13I14H loose their conventional 
interpretability for (3 < da in d dimensions. The value, 
Th = 1/ (da) is a limiting temperature for positive val- 
ues of the parameter a, at which the physically relevant 
quasi-energy per particle diverges, so there is no use of 
further heating at this temperature. More and more en- 
ergy given to the system would raise the temperature less 
and less. For a < 0, i.e. for attractive correction, the en- 
ergy per particle is limited by (1 + 1/2 + 1/3)E C from 
above, but any temperature may occur (cf. Fig^]). 

The ./V-particle energy distribution can be obtained 
also exactly in this case. The Fourier-integral 10 has 
TV-fold poles at the values Sk = — iBk for each direc- 
tional degree of freedom k = 1, . . . , d. The evaluation of 
such integrals is somewhat involved in a general number 
of dimensions, so we restrict our further analysis to the 
cases d = 1 and d = 3. For a one-dimensional Tsallis 
distribution we obtain 



g N +i(X) = AEB\ 1 J exp(X/E c 



X/T). (15) 



In this case a Hagedorn exponential emerges with the 
limiting temperature Th = E c . Considering on the other 
hand F N+1 (E) = g N+1 (X)X'(E), i.e. the N + 1-particlc 
distribution of the (not conserved) naive energy expres- 
sion only, the exponetially growing factor does not occur. 



It is due to X'(E) 



-aX 



for the Tsallis distribution: 



, , Bf +1 (ln(l 
F N+1 (E) = AE^- 



aE)) 



N 



N Nl 



-(l + aE)- 1 /^. 



(16) 
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FIG. 1: The temperature, T, as a function of the quasi- 
energy per particle, (X) /N , in units of the Tsallis energy scale 
E c = 1/a. The linear relation E/N = 3T known for massless 
particles is realized for a = 0. For repulsive corrections, a > 0, 
a limiting temperature occurs, Th = l/3a = E c /3. 



However, still j3 > a or T < Th must be satisfied, and 
for positive a values increasing quasi-energy per particle 
does not raise the temperature above Th- 

The corresponding expression in d = 3 dimensions is 
more involved. There occur three different exponential 
factors with lowest limiting temperature Th = E c /3. 
The general dependence on (5 can be factorized out by 
shifting the integration variable s to s + i[3: 

yN+1 

9N+l{X) = {^Zmr^ '^^^ 

with $>n(X) depending on the dynamical input parame- 
ters only, but not on the temperature T = 1//3: 



N+l 



sA 



nv(< 



(18) 



\k=l 



Factorizing out the ideal thermal factor e~" x /Zn, the 
rest can be regarded as the mass spectrum of the N- 
parton system at X = m and AE = Am: 



p N (m) 



V 



N 



$jv_i(m). 



(19) 



The functions $jy(X) obey the recursion rule 



$ N (X)= [ dt$ N ^(t)<f> (X-t) (20) 
Jo 

with the starting point of the recursion, 

§ (X) = (e 3aX - 2e 2aX + e aX ) . (21) 

Particular important cases are N = (partons), N = 1 
(mesons or diquarks) and N — 2 (direct baryons). Be- 



FIG. 2: The integrated, non-degenerated mass spectrum of 
mesonic and baryonic resonances from the 2004 compilation 
by the Particle Data Group |30lj is compared with the respec- 
tive cumulative numbers of two and three-parton clusters. 



sides the already known ^q(X) (cq. I21fl we obtain 

(aX - 3)e 3aX + AaXe 2aX + (aX + 3)e aX 



and 



4a 5 



(22) 



16a s 



[((ol) 2 - 9aX + 24) e 3a 



x 



2 i fi\ „2aA 



-8((aX) 2 + 6 
+ {{aXf + 9aX + 24) e aX ] (23) 



For any TV the result contains three exponentials giving 
rise to a lowest limiting temperature of Th = E c /3 = 
l/3a for positive values of the parameter a. 

Fig[21 shows the non-degenerate, integrated mass spec- 
tra. Values published on the Particle Data Group home- 
page are summed up in mass histograms. The respective 
numbers are fitted as Nm = l+A/i((m-mM)/3T/f) and 
N B = 1 + A 2 f 2 ((m - m B )/3T H ) with A = V/a 3 ir 2 and 
the integral functions f n (x) = <b n {t)dt. The fits to the 
data are most sensitive to the value of Th, which however 
may be compensated by changing the assumed volume, 
V. Keeping wlm = 0.14 GeV and tub = 0.94 GeV, only 
A and Th are varied. Our best fit gives a relatively high 
value, T H = 0.35 GeV (meaning E c = 1.05 GeV) and 
a volume of V — 261 /to 3 (a sphere with a radius of 4 
fm, or a box with a length of 6.4 fm). Above the masses 
where the data seem to deviate from the fast growing 
part, the fit cannot be followed any more. According to 
ref.[2j, newest data raise the experimental curve higher. 
Our idea, different from both the string and bag model 
consideration, seem to agree with the difference between 
the meson and baryon mass spectra, as well as with a 
polynomial upcurving of the baryon spectrum. 

In conclusion wc pointed out that Tsallis distributed 
massless partons can be combined to an effective mesonic 
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and baryonic mass spectrum by considering the con- 
served quasi-energy as the hadron energy. Besides an 
ideal thermal factor, e _/3X /Zjf, a further energy depen- 
dent factor results from the folding of parton distribu- 
tions. It can be regarded as a thermal (pre)hadron 
mass spectrum emerging from a statistical hadroniza- 
tion picture. The prediction of this folding, while having 
two parameters (a volume and the Hagcdorn tempera- 
ture) , gives an acceptable qualitative agreement with the 
known hadronic mass spectrum. In this picture a natural 
difference emerges between mesonic and baryonic reso- 
nances due to their different foldness by parton coales- 
cence. The characteristic temperature, Tu = l/3a ~ 350 
MeV is a limiting temperature: one cannot increase the 
temperature above this value, not even with an infi- 
nite amount of energy. The parameter E c = 1/a = 
3Th ~ 1-05 GeV provides at the same time the scale 
where the power-law tail of individual p^-spectra starts 
to dominate the exponential part, and it is intimately 
related to the typical pair-interaction energy, due to 
h(Ei,E2) — Ei — E-2 = EiE2/E c . 
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